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Abstract  for  Poster 

for  the  4th  Annual  University  of  California 
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26-28  March  1990 

Avalanche  Dynamics  and  Self-Organization  in  a  Magnetic  Froth 
K.L.  BABCOCK*  and  R.M.  WESTERVELT,  Harvard  University. 

Uniaxial  magnetic  garnet  films  support  an  experimentally  accessible  system  of  two- 
dimensional  domain  patterns  that  evolve  in  response  to  external  magnetic  fields.  This  poster 
will  describe  observations  of  the  coarsening  of  disordered  cellular  domain  patterns  in  an 
increasing  bias  field.1  The  patterns  coarsen  via  the  shrinking  and  elimination  of  cells  with 
fewer  than  six  sides,  a  process  topologically  similar  to  that  observed  in  soap  froths  and  the 
annealing  of  polycrystalline  materials.  Novel  to  the  domain  patterns  are  structures  containing 
trapped  magnetic  bubbles  which  act  as  stable  five-fold  vertices.  These  "bubble  traps"  resist 
elimination,  and  so  provide  a  barrier  to  the  coarsening  that  allows  the  patterns  to  settle  into 
stationary,  metastable  states.  Bias  increments  destabilize  these  states,  leading  to  avalanche¬ 
like  reorganizations  in  which  cells  are  destroyed  by  a  domino  effect.  The  avalanche 
dynamics  are  self-organizing:  for  significant  ranges  of  bias  field  and  cell  density,  the 
stationary  patterns  are  barely  stable  with  respect  to  small  bias  increments.  The  avalanche 
sizes  and  lifetimes  span  more  than  two  decades  and  have  distributions  that  are  best  described 
as  power  laws  over  a  limited  range.  These  features  suggest  that  this  is  an  experimental 
manifestation  of  self-organized  criticality.2  The  self-organizing  dynamics  also  lead  to  scaling 
in  the  structure  of  the  patterns:  the  fraction  of  cells  that  are  bubble  traps  remains  nearly 
constant  over  a  two-decade  drop  in  cell  density.3 

*  Currently  at  the  University  of  California  at  Santa  Barbara,  Dept,  of  Physics. 

1  K.L.  Babcock  and  R.M.  Westervelt,  Phys.  Rev.  A  40,  2022  (1989). 

2P.  Bak,  C.  Tang,  and  K.  Wiesenfeld,  Phys.  Rev.  Lett.  59,  381  (1987). 

3  K.L.  Babcock  and  R.M.  Westervelt,  Phys.  Rev.  A  41,  Feb.  15  (1990). 


THE  SHAPES  OF  RANDOM  WALKS  WITH  FIXED  END-TO-END  DISTANCE 


Arezki  Beldjenna  and  Joseph  Rudnick 

Department  of  Physics 
University  of  California  at  Los  Angeles 
Los  Angeles,  CA  90024 

and 

George  Gaspari 

Department  of  Physics 
University  of  California  at  Santa  Cruz 
Santa  Cruz,  CA  95064 

ABSTRACT 

Anisotropies  of  random  walks  embedded  in  d  spatial  dimensions  and  with  fixed  end-to-end 
distance  are  discussed.  As  a  representation  of  polymers,  these  walks  correspond  to  long  chain 
molecules  whose  heads  and  tails  are  fixed  in  space.  An  exact  analytical  expression  for  the 
asphencity,  valid  in  an  arbitrary  spatial  dimensionality,  is  presented.  We  also  present  expressions 
for  the  average  principal  radii  of  gyration  to  first  order  in  lid.  The  results  reveal  two  regimes  that 
depend  on  the  fixed  length.  They  recover  the  results  for  both  i  nrestricted  open  and  closed  random 
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Determination  of  Excess  Properties  and  infinite  Dilution 
Quantities  fron  Computer  Simulation  via  Single 
Charging-Integral  Approach 

by 

Ariel  A.  Chialvo  and  Pablo  G.  Debenedetti 
Department  of  Chemical  Engineering 
Princeton  University 
Princeton,  NJ  08544-5263 

Rigorous  statistical  mechanical  expressions  are  derived  to  evaluate  activity 
cceffic.ents,  partial  molar  volumes  and  partial  molar  enthalpies  at  infinite 
dilution  via  NPT  molecular-based  computer  simulation.  These  equations  are 
based  on  the  recently  proposed  single  charging-integral  method  for  the 
determination  of  excess  properties  in  multicomponent  model  systems  [1]. 

The  distinguishing  feature  of  the  new  method,  in  contrast  to  conventional 
ones,  is  that  there  is  no  need  for  simulations  at  infinite  dilution  or  for  unreliable 
extrapolation  in  order  to  obtain  infinite  dilution  quantities. 

Since  the  single  charging-integral  method  is  based  on  the  coupling 
parameter  concept,  the  approach  avoids  all  of  the  problems  commonly 
associated  with  particle  insertion  methods. 

Simulation  results  for  infinite  dilution  activity  coefficients,  Henry  constants 
and  pure  component  fugacities  of  a  binary  Lennard-Jones  mixture  are 
presented  and  compared  with  theoretical  results  from  two  perturbation  schemes 
[2,3] 


[1]  Chialvo,  A. A.,  J.  Chem.  Phys.,  In  press  (1990) 

[2]  Shukla,  K.P.,  Molec.  Phys.,  £2,  1143  (1987) 

[3]  Fischer,  J.  and  S.  Lago,  Chem.  Phys.,  Zfl.  5750  (1983) 


Shear  Viscosity  Calculation  via  Equilibrium  Molecular 

Dynamics: 

Einstenian  vs.  Green-Kubo  Formalism 

by 

Ariel  A.  Chiaivo  and  Pablo  G.  Debenedetti 
Department  of  Chemical  Engineering 
Princeton  University 
Princeton,  NJ  08544-5263 

Twenty  years  after  the  first  report  on  shear  viscosity  calculation  via  Moiecular 
Dynamics  [1],  there  is  still  controversy  over  the  proper  evaluation  of  this  property 
[2],  An  intriguing  fact  is  that  despite  the  numerical  advantages  of  the  Einstenian- 
like  (i.e.,  squared  displacement)  formalism,  most  (if  not  all)  work  has  been  done 
through  the  application  of  the  Green-Kubo  approach.  Although  the  theoretical 
equivalence  between  both  formalisms  was  demonstrated  by  Helfand  [3],  their 
relative  numerical  and  computational  merits  are  quite  different.  Furthemore,  the 
mean  squared  displacement  expression  can  be  written  in  two  theoretically 
equivalent  but  numerically  very  different  forms:  one  given  by  Helfand  [3],  and  a 
much  less  frequently  used  form  given  by  McQuarrie  [4], 

In  this  work,  we  compare  the  numerical  and  computational  implications 
(accuracy,  memory  requirement,  sample  size  and  cutoff-dependence)  of  the  two 
Einstenian-like  equations  and  of  different  expressions  based  on  the  Green- 
Kubo  formalism.  In  contrast  to  Hoheisel  and  Vogelsang’s  conclusions  [2],  we 
find  that  the  Einstenian  formalism  gives  excellent  results  and  is  less 
computationally  demanding  than  its  Green-Kubo  counterpart.  In  addition,  we 
find  that  the  expression  given  by  McQuarrie  [4]  is  numerically  more  convenient 
than  Helfand's  expression  [3],  since  it  requires  a  smaller  number  of 
"experiments"  (uncorrelated  time  origins).  All  of  our  conclusions  are  validated 
by  extensive  comparisons  with  literature  values  for  shear  viscosity  of  Lennard- 
Jones  systems. 

[1]  Alder,  B.J.  D.M.  Gass,  and  T.E.  Wainwright,  J.  Chem.  Phys.,52,  3813 
(1970) 

[2]  Hoheisel,  C.  ,  and  R.  Vogelsand,  Comp.  Phys.  Rep.,  £,  1  (1988) 

[3]  Helfand,  E.,  Phys.  Rev.,  Hi,  1  (1960) 

[4J  McQuarrie,  D.  ,  "Statistical  Mechanics",  Harper&Row,  New  York  (1976) 
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Molecular  Bonding  and  Packing  of  Sulfur 

L.  Rene  Corrales  and  Peter  J.  Rossky 


A  novel  format  for  modeling  the  formation  of  bonds  and  packing  arrangements  in 
covalently  bonded  materials  is  introduced  and  applied  to  sulfur.  The  method  takes  into 
consideration  the  implicit  contributions  from  chemical  bonding  and  conventional  packing 
interactions  using  ideas  borrowed  from  both  semi-empirical  quantum  chemistry  and 
molecular  mechanics.  This  method  has  the  computational  advantage  of  being  efficient, 
inexpensive  and  broadly  applicable.  Results  for  sulfur  are  presented  and  contrasted  with 
experimental  data  and  other  theoretical  methods. 


wader  f\"<=  ne  ^  e-'m  i  red  .  w  i  “h  m  th  --  nc>c  -  .  tv.  rhe  c; -  r  :  u  ~  -  ■.  r  .  ■ 

.■if*  J3fjrj  Ti'iC-  7.  ■.  7  -7  r  ,  :  ■  ?  i.  f  ! ‘iM  *-  i;  1  .  ■-  -  ;  _•  ;  ■: :  .  ;■ 

tr.vns  i  t  ion :  the  nematic  i.auid  :  s  coos  ice-  -rz  a-r  \  o  :  v  t 

in  which  each  component  reprsiert  •  a  sc  c :  i  *'  :c  .nr-  leculi, 
o  r  l  en  tali  o  n  .  We  found  that  the  f  ir?t  order  hr  a  ns  i  t ;  on  can  be 
recovered  through  a  trier  jticai  or  a  critical  e  id  point.  b/ 
applying  an  external  field  to  the  neta'cc  phase. 

h  mean  field  lattice  model  is  used  to  predict  the  phase 
diagrams  of  a  system  that  forms  micro emu  1 s ions  in  different 
physical  situations.  The  calculated  phase  progressions  and  toe 
mterfacial  proper  •;  ies  agree  qua!  i  datively  with  etpsr  -.mental 
results.  We  conclude  that  the  ultra-low  inter f ac lal  tensions  we 
get  are  due  to  the  presence  of  a  tricritical  point  near  the  three 
phase  line.  We  also  discuss  the  phase  diagrams  for  systems  that 
show  ordered  phases. 
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Ve  analice  some  hvpotbetical  conditions  ir.  which  rancora 
fluctuations  of  ion  concentration  in  restricted  regions  sight 
lead  to  transient  ion  flow  through  a  cell  or  mitochondrial 
membrane,  and  to  eventually  produce  ion  concentrations  gra¬ 
dients  across  these  structures.  We  observe  from  figures  that 
the  larger  the  required  relative  fluctuation  £  ,  the  smaller 
is  the  corresponding  volume  enclosing  it.  However,  there  is  a 
broad  range  of  concentrations  for  which  significant  fluctua¬ 
tions  of  c  cake  place  in  enclosing  volumes  having  diameters 
corcoarabie  to  the  membrane  thickness.  In  these  cases,  the 
dissipation  of  a  fluctuation  close  to  the  membrane  could  re¬ 
sult  in  a  transient  net  flow  of  the  corresponding  ion  through 
the  channel.  In  the  case  of  protons  the  fluctuation  in  con¬ 
centration  should  have  considerable  significance,  due  to  the 
very  low  basal  H.+-  concentration,  tor  example,  in  a  volume  of 
one  cubic  micrometer  there  are  about  60  free  protons  at 
pH  ?  and  this  number  should  fluctuate  between  52  and  66.  what 
should  produce  significant  changes  in  the  local  pH. 

(Partially  supported  by  CNPq-Bratil  grant  200-63/89-6/EP} . 


,  c  •• 

n  uw.tuatio»»  tut  (m  ’Tt  >i  i  ittyt  rum 

trritAt.  n*  co<>CttitBAT  luhs  <c>  ahu  i*u  or  ■ioiocicaw 

rio»oi«ic  riurtUAi  Hint . 


MIU  UUtlKt  Or  TAkliniJ  flivt>l*r®  |«)  n 

IjjUi;.  {.)  10"»C  I «  U«."Tl«AT  I0t<  IM  UJilMi  \ 


BINARY  MIXTURES  IN  SILICA  GELS 


B.J.  Frisken,  F.  Ferri  and  D.S.  Cannell 
Department  of  Physics,  UCSB 


We  have  studied  the  pretransitional  behavior  of  the  binary  mixture  of  lutidine 
and  water  in  silica  gels.  The  gels  used  were  made  by  polymerization  of  TMOS 
(tetramethylorthosilicate)  in  water.  Static  light  scattering  studies  show  that  these 
gels  are  fractal  on  length  scales  ranging  up  to  a  crossover  length  which  can  be 
varied  from  200-4000  A  by  gelation  conditions.  At  length  scales  greater  than  the 
crossover  length,  the  gels  act  as  a  randomly  distributed  collection  of  fractal 
objects.  Samples  were  prepared  so  that  the  lutidine  volume  fraction  in  the  gels 
was  within  1%  of  that  expected  in  a  sample  at  bulk  critical  concentration.  The 
behavior  of  the  mixture  in  the  gel  in  the  one  phase  region  was  then  studied  by 
ftatic  light  scattering.  It  was  found  that  the  distribution  of  the  intensity  of 
scattered  light  with  scattering  vector  q  resembled  the  form  of  that  of  the  bare  gel, 
with  an  amplitude  at  q=0  which  showed  power  law  divergence  as  the  transition 
temperature  of  the  bulk  mixture  was  approached.  Possible  explanations  for  this 
power  law  behavior  will  be  discussed. 


Fluctuating  Hydrodynamics  and  the  Principal  Oscillation  Pattern  analysis 


Alejandro  Garcia 

Dept.  Physics,  San  Jose  State  Univ.,  San  Jose  CA  95192-0106 

Cecile  Penland 

CIRES,  Univ.  of  Colorado,  Boulder  CO  80309-0216 

Principal  Oscillation  Pattern  (POP)  analysis  was  recently  introduced  in 
climatology  to  analyze  timeseries  xgt)  produced  by  systems  whose 
dynamics  is  described  by  the  stochastic  differential  equation,  xj  =  Bij  xj  + 
c,i.  The  matrix  B  gives  the  deterministic  feedback  and  q  is  a  white  noise 
with  covariance  <qi(t)qj(t’)>  =  Qij  S(t-t').  The  POP  method  was  applied  to 
data  from  a  Direct  Simulation  Monte  Carlo  program.  The  system  was  a 
dilute  gas  with  50,000  particles  in  a  Rayleigh-Benard  configuration.  The 
POP  analysis  correctly  obtained  the  feedback  matrix,  B,  as  given  by  the 
Navier-Stokes  equations  and  the  covariance  matrix,  Q,  as  given  by  Landau- 
Lifshitz  fluctuating  hydrodynamics.  We  discuss  how  the  method  may  be 
used  to  empirically  compute  amplitude  equation  coefficients  from  data. 


series  investigations  of  frustrated  square-lattice  antiferromagnets 


The  zero-cemperature  properties  of  the  square-lattice,  S= 1/2 
Heisenberg  model  with  antiferromagnetic  first  and  second  (diagonal) 
neighbor  interactions  have  proven  to  be  remarkably  elusive.  For 
J.SJdJ.,  ,  there  are  at  least  three  proposals  for  the  nature  of  the 
ground  state  (disordered,  helical,  spontaneously  dimerized),  and 
rinite-size  studies  have  vielded  conflicting  conclusions.  We  report 
on  series  expansions  in  coupling  constants  about  a  variety  of 
trivial  (explicitly  dimerized)  Hamiltonians.  The  calculations 
suggest  that  in  the  disputed  regime,  the  model  is  spontaneously 
dimerized,  buc  only  weakly  so.  Furthermore,  there  appear  to 
be  several  relatively  large  correlation  lengths. 


SURFACE  DIFFUSION  AND  FLUCTUATIONS  OF  GROWING  INTERFACES 

Leonardo  Gol ubovi c  and  Robijn  Bruinsma 
Physics  Department  and  Solid  State  Science  Center 
University  of  Cal i f orn i a , Los  Angeles, CA  90024 


Abstract 

A  model  is  proposed  for  the  growth  of  realistic  interfaces 
relaxing  dominantly  by  surface  diffusion.  We  show  that  the  thermal 
noise  can  produce  large  scale  orientational  fluctuations  in  the  active 
zone  of  the  growth.  We  relate  them  to  the  experimentally  observed 
undersurface  networks  of  voids  which  we  characterize  by  a  structural 
lengthscale  of  the  order  of  the  orientational  persistence  length. 


Roughening  transition  of  the  isotropic  frustrated  Ising  model 


B.  Kahng  and  K.  A.  Dawson 

Department  of  Chemistry,  P.  0.  Box  224,  University  of  California 

Berkeley.  CA  94720 

Abstract 

We  study  the  roughening  transition  behavior  with  the  isotropic  competing  Ising  Hamil¬ 
tonian  1  ,  H  =  -J  Ln  *i*}  -  2M  Tiinn  s tSj  -  M  Vnrin  s,s;.  where  mi,  dnn.  nnn  mean 

the  nearest  neighbor,  the  diagonal  neighbor,  the  linearly  next  nearest  neighbor  respec¬ 
tively.  and  J  >  0  and  —  x:  <  \I  <  .  Low  temperature  expansions  are  performed  in  the 

three  dimensional  cubic  lattice  for  moments  of  the  interface  width  and  for  the  interface 
density.  The  results  are  compared  to  a  Sine-Gordon  renormalization  group  calculation  of 
the  analogous  solid-on-solid  model.  By  comparison  with  Monte  Carlo  simulation  results 
U  it  is  likely  that  the  roughening  transition  line  in  the  (j  =  J/kT.m  =  M/kTi  plane 
passes  through  the  Lifshitz  point,  and  is  very  close  to  part  of  the  phase  boundary  between 
paramagnetic  and  modulated  phases.  We  discuss  the  implications  for  the  phase-structure 
and  the  interfacial  profile  in  the  vicinity  of  the  Lifshitz  point. 
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About  negative  moments  of  current  distribution 
in  random  resistor  networks 

Byungnam  Kahng 

Department  of  Chemistry 
University  of  California,  Berkeley,  CA  947201 
Center  for  Polymer  Studies  and  Department  of  Physics 
Boston  University.  Boston  \IA  02215 


Abstract 

We  investigate  the  size  dependent  behavior  for  negative  moments  of  current  distribu¬ 
tion  in  random  resistor  networks.  With  the  extreme  value  statistics  idea,  we  derive  the 
explicit  form  of  the  current  distribution  in  small  current  regime.  The  distribution  is  found 
to  be  the  Weibull  distribution  n(x)  ~  exp(— xm),  with  m  >  1  and  x  =  —  logz/log£.  £ 
is  the  system  size.  Using  this  distribution  function,  we  show  breakdown  of  multifractal 
behavior  in  negative  moments  regime  and  non -occurrence  of  thermodynamic  phase  tran¬ 
sition. 

P.A.C.S.  numbers:  64.60. Ak. 05.50. -4-q, 05.60. -t-w 
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MUSIC  COMPOSITION  AND  MODELS  OF  CORTEX 
USING  PHYSICAL  SYSTEM  ANALOGIES 

Xiaodan  Leng  and  Gordon  L.  Shaw,  Center  for  the  Neurobiology  of  Learning  and 
Memory  and  Physics  Department  University  of  California,  lrvine,  Irvine,  CA  91717. 
Enc  Wright  Irvine  Conservatory  of  Music,  711  Farad  St.  Costa  Mesa,  CA  91627 

Based  on  the  Fisher’s  ANNNI  model  of  physical  spin  system  and  guided  by 
the  Mountcastle  cortical  principle  of  columnar  organization  we  had 
developed  the  trion  model  [  G.  Shaw,  D.  Silverman  and  J.  Pearson,  PNAS 
82,2364  (1985)]  for  the  firing  activity  of  groups  of  neurons;  an  essential 
feature  of  the  trion  network  is  that  it  is  highly  structured  in  time  and  in 
spatial  connections.  The  trion  model  has  been  shown  to  be  useful  for 
understanding  aspects  of  memory  and  pattern  recognition.  To  study  the 
properties  of  interconnected  networks  of  the  trion  model  we  incorporated 
the  MILA  model  [J.  Patera,  G.  Shaw,  R.  Slansky,  and  X.  Leng,  Phys.  Rev. 
A40,  1073  (1989)  whose  underlying  structure  is  the  weight  systems  of 
zero-level  representations  of  affine  Kac-Moody  algebras]  to  modulate  the 
firing  activity  of  interacting  cortical  columns.  Turning  now  to  higher 
cognitive  processes,  we  map  spatio-temporal  firing  pattern  evolutions  of 
the  model  onto  pitches  and  instruments  to  produce  music.  The 
preliminary  results  are  striking:  1)  The  trion  music  is  very  pleasant  to 

listen  to:  I'l  On*  c«f  of  connections  among  the  trions  produces  music 

sounding  like  a  waltz,  another  set  gives  a  minuet;  3)  Different  mappings  of 


Monte  Carlo  Calculations  of  the  Correlation  Functions  for 
Heisenberg  Spin  Chains  at  T=0 


Shoudan  Liang  ,  Department  of  Physics  and  Material  Research 
Laboratory,  University  of  Illinois  at  Urbana-Champaign,  1110 
West  Green  Street,  Urbana,  Illinois  61801. 

Abstract:  We  computed  the  ground  state  spin-spin  correlation  function  for 
the  one  dimensional  antiferromagnetic  Heisenberg  model  of  spin  1/2,  1  and 
3/2.  The  spin-1  chain  has  a  disordered  ground  state  with  a  correlation 
length  of  6.2  lattice  spacing  whereas  the  spin- 1/2  and  the  spin-3/2  chains 
both  have  (ln(r))a / r  decaying  correlation  functions.  The  logarithmic 
corrections  are  different  for  the  spin- 1/2  and  the  spin-3/2  chains.  Finite-size 
scaling  is  used  to  analyze  the  data  for  chains  of  up  to  128  spins. 


Existence  of  N^ei  Order  at  T=0  in  the  spin-1/2 
antiferromagnetic  Heisenberg  Model  on  a  Square  Lattice 

Shoudan  Liang  .  Department  of  Physics  and  Material  Research 
Laboratory,  University  of  Illinois  at  Urbana-Champaign,  1110  West 
Green  Street,  Urbana,  Illinois  61801. 

Abstract-.  We  have  computed  spin-spin  correlation  function  for  the 
ground  state  of  the  spin- 1/2  antiferromagnetic  Heisenberg  model  on 
a  32x32  square  lattice  using  a  projection  Monte  Carlo  method. 

Kennedy,  Lieb  and  Shastry  have  shown  that  if  the  ground  state  is  not 
N'<iel  ordered,  the  two  point  function  must  satisfy  an  upper  bound. 
Our  correlation  function  at  short  distance  (less  than  4.5)  violates  this 
bound  which  implies  that  the  ground  state  must  have  long  ranged 
N6el  order.  We  have  also  estimated  the  staggered  magnetization  to 
be  rrL=0.304±0.04  in  the  thermodynamic  limit. 


ABSTRACT  FOR  SUBMISSION  TO  UC  STATISTICAL  MECHANICS  CONFER¬ 
ENCE 


THE  DIFFUSION  COEFFICIENT  OF  A  REPTATING  POLYMER.  J.M.  Deutsch  and 
T.L.  Madden.  University  of  California  at  Santa  Cruz,  Santa  Cruz,  CA,  95064. 


The  dynamics  of  a  polymer  in  a  network  of  entanglements  is  studied.  The  viscosity 
and  diffusion  coefficient  were  studied  in  the  repton  model  for  up  to  75  links.  First  we 
derive  the  correct  procedure  for  calculating  the  three  dimensional  diffusion  coefficient 
in  this  model.  We  find  that  the  diffusion  coefficient  scales  as  £--~5±o.i.  This  result  is 
in  marked  contrast  to  the  asymptotic  L~ 2  prediction  of  reptation.  We  see  that  the 
asymptotic  behavior  of  the  diffusion  coefficient  must  have  corrections  of  the  order  of 
L~v 2.  This  implies  that  fluctuations  in  the  tube  length  cannot  provide  an  adequate 
explanation  for  the  experimental  data.  The  viscosity  is  found  to  scale  as  T3-49±013t  in 
accord  with  experimental  data  and  other  numerical  results. 


Abstract  for  the  UC  Conference  on  Statistical  Mechanics 
Santa  Barbara,  Ca.  March  1990 

PHASE  TRANSITIONS  IN  THE  TRION  MODEL 
OF  CORTICAL  ORGANIZATION 


John  McGrann,  Gordon  L.  Shaw,  Dennis  J.  Silverman 
Department  of  Physics,  University  of  California,  Irvine 
Irvine,  Ca.  92717 

John  C.  Pearson 

David  Sarnoff  Research  Center,  Princeton,  NJ  08543 


The  trion  model  [Shaw,  Silverman,  Pearson,  PNA3 
82,2364,(1985;],  based  on  an  analogy  with  Fisher's  ANNXI  physical 
spin  system  [Fisher, Selke,  Phys.  Rev.  Lett. ,44. (1980)],  has  been 
motivated  by  such  cortical  principles  as  Mountcastle's  columnar 
organization.  The  trion  represents  a  localized  group  of  neurons  with 
three  levels  of  firing  activity  which  has  spatial -temporal  connections 
that  are  highly  structured.  Networks  of  a  small  number  of  trions 
(with  structured  interactions  and  firing  dependent  on  activity  at  two 
previous  time  steps)  support  a  repertoire  of  hundreds  of  quasi- 
stable,  periodic  faring  patterns,  M?s  (which  can  be  learned).  We 
report  striking  phenomena  with  variation  in  "temperature'’  T:  There 
exist  a  series  of  "phase  transitions"  at  precise  values  T(n)  giving  new 
repertoires  of  MPs,  and  the  average  time  for  any  initial  firing 
configuration  to  project  onto  an  MP  shows  a  quite  sharp  change  at 
each  T(n).  When  the  connections 
are  made  random  these  phenomena  disappear. 


Symmetry  Breaking  in  a  Simple  Electron-Phonon  Model  with  Electron 

Correlations 

Qing  Zhang*  and  Robert  J.  Silbey 
Department  of  Chemistry 
Massachusetts  Institute  of  Technology 
Cambridge  MA  02139 

Abstract 

We  consider  a  simple  model  which  contains  both  electron  phonon  coupling, 
electron  transfer,  and  electron  correlation,  which  can  be  a  prototype  of  a  mixed 
valence  system  with  3  nearby  electronic  states.  We  treated  the  adiabatic 
solution,  simple  corrections  to  that,  and  the  effect  of  symmetry  breaking  terms 
on  the  solutions.  We  find  that  the  adiabatic  solutions  in  the  symmetric  case  can 
show  a  (false)  symmetry  breaking  which  is  removed  in  the  nonadiabatic 
calculation.  These  solutions  are  extremely  unstable  with  respect  to  the  addition 
of  a  nonsymmetric  term  in  the  Hamiltonian. 


*  Current  address:  Department  of  Chemistry,  UC  Berkeley 

Berkeley,  CA  94720 


ABSTRACT 


?o$  /e-K' 

Title  of  Diiatrtacion:  EQUILIBRIUM  POLYMERIZATION  AS  A 

PHASE  TRANSITION 

Jaime  Ruiz  Garcia,  Doctor  of  Philosophy,  1989 

Dissertation  directed  by:  Sandra  C.  Greer 

Professor 

Department  of  Chemistry 

Recently,  Wheeler  et  al.  have  proposed  a  theory  to 
describe  equilibrium  polymerization  reactions.  The  theory 
is  based  on  the  formal  limit  n  ->  0  of  the  n-vector  model 
of  magnetism.  Wheeler  et  al .  have  predicted  that  the 
equilibrium  polymerization  reaction  can  be  treated  as  a 
phase  transition,  analogous  to  the  phase  transitions 
observed  in  magnetic  systems. 

Several  predictions  of  Wheeler  et  al .  have  been 
experimentally  tested  in  sulfur  and  sulfur  solutions, 
e  g. ,  the  phase  diagrams,  the  divergence  of  the  heat 
capacity  at  the  polymerization  transition,  etc..  There 
are  some  discrepancies  between  the  theory  and  the  ex¬ 
perimental  data  available  in  the  literature  for  the 
sulfur  systems  . 


evidence  of  Che  presence  of  a  tricritical  point.  However, 
some  discrepancies  still  exist.  We  have  not  a  good 
explanation  for  these  differences;  thus  more  work  is  still 


needed . 
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HEXAGONAL  AND  DISORDERED  PHASES  OF  SEMIFLEXIBLE  CHAINS. 
I.  V.  Selinger  and  R.  F.  Bruinsraa,  UCLA.* — We  consider  the  statistical 
mechanics  of  a  system  of  semif lexible  chains,  which  can  represent 
polymer  liquid  crystals,  long-chain  biomolecules,  stiff  worm-like 
micelles,  or  columns  of  discotic  liquid  crystals.  In  the  orderered 
phase,  these  systems  can  form  a  hexagonal  lattice;  in  the  disordered 
phases,  they  can  become  nematic  or  isotropic  fluids.  We  calculate  the 
structure  factor  and  the  orientational  correlation  tunction  in  tne 
hexagonal  and  nematic  phases  using  continuum  elastic  theory.  In  the 
hexagonal  phase,  the  structure  factor  exhibits  Bragg  peaks  surrounded  by 
diffuse  thermal  scattering.  The  diffuse  line  shape  is  related  to  the 
elastic  constants  of  the  lattice.  In  the  nematic  phase,  the  structure 
factor  exhibits  a  similar  line  shape  as  the  wavevector  q— >0.  In 
general,  the  elastic  constants  depend  on  the  type  of  interaction  between 
the  chains.  We  derive  these  constants  for  a  pure  excluded-volume 
repulsion.  We  then  use  a  variational  theory  to  model  the  first-order 
transitions  among  the  hexagonal,  nematic,  and  isotropic  phases.  This 
theory  leads  to  a  phase  diagram  in  terms  of  pressure  (or  concentration) 
and  persistence  length. 

*Supported  by  the  UCLA  Solid  State  Science  Center. 
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Abstract 

A  simple  novel  scaling  argument  for  diffusion- limited  reactions  for  both  A  +  A  — >  0  and 
A  +  b  — >  0  is  developed.  Asymptotic  behaviors  of  species  densities  without  source  input  and 
reaction  orders  with  source  input  are  discussed.  It  is  shown  how  the  structures  of  spatial  dis¬ 
tributions  of  particles  affect  reaction  rates.  As  a  result,  the  asymptotic  behavior  of  the  density 
for  A  +  B  — »  0  without  source  input  is  found  to  depend  on  the  fractal  dimension,  a  depen¬ 
dence  that  does  not  seem  to  have  been  noted  before. 
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Thermodynamic  Equations  of  State  for  the  Kinetic 
and  Potential  Energies  of  a  System  of  Charged  Particles 


DUANE  P.  SIEMENS 
Department  of  Physics 
University  of  California 
Davis,  C.4  95616 


ABSTRACT 

In  this  paper  we  consider  a  system  of  N  charged  particles  confined  to  a  volume  V . 
We  prove  the  following  relationship  between  the  thermodynamic  averages  of  the  kinetic 
energy.  <  T  >,  and  potential  energy,  <  V  >, 

2  <  T  >  -  <  V  >  -3PV  =  0, 

where  P  is  the  pressure.  Using  this  relationship  we  relate  the  average  kinetic  and  potential 
energies  to  measurable  thermodynamic  quantities.  We  also  derive  an  analogous  realtion- 
ship  for  a  given  species  of  particle  in  a  multi-species  system.  We  further  consider  the  effect 
of  a  constant  external  magnetic  field. 


PHASE  TRANSITION  BEHAVIOR  IN  A 
NEGATIVE  TEMPERATURE 
GUIDING  CENTER  PLASMA 
OR  VORTEX  GAS1 

Ralph  A.  Smirh 

Physics  Dep’t.  B-019.  LCSD.  La  Jolla.  CA  92093 

The  thermal  equilibria  of  a  two-dimensional,  single-species,  guiding- 
center  plasma  or  gas  of  parallel  line  vortices,  bounded  by  a  circular 
equipotential  surface,  are  considered  in  the  microcanomcal  ensem¬ 
ble.  The  system  is  a  model  for  a  magnetically  confined  nonneutral 
plasma,  or  an  extended  vortex  column  in  an  ideal  fluid.  At  sufficiently 
large  energies  i  for  specified  canonical  angular  momentum),  the  statis¬ 
tically  predominant  configurations  are  displaced  macroscopically  from 
the  axis.  The  transition  between  symmetric  and  asymmetric  equilib¬ 
ria  resembles  a  second-order  phase  transition,  and  occurs  at  formally 
negative  temperatures.  It  is  related  to  a  bifurcation  in  the  mean-field 
description.  An  analysis  of  fluctuations  suggests  that  the  transition 
is  of  the  Gaussian  type,  in  the  limit  of  large  particle  number.  The 
theory  is  compared  with  microcanonical  Monte  Carlo  simulations. 


:  Work  supported  by  ()  N  R  N000 11-82- K-06'21  ami  N  S  F  PH Y87-0t>358;  Computa¬ 
tions  done  at  die  San  Diego  Supercomputer  (.’enter 


UNIVERSALITY  IN  SOLUTIONS  OF  LINEAR  POLYMERS’ 


Lipiin  Sung  and  David  S.  Canned 

Department  of  Physics,  University  of  California,  Santa  Barbara,  CA  93106 


The  universal  properties  of  linear  polymer  solutions  have  been  intensively  investigated 
both  experimentally  and  theoretically.  For  appropriate  solutions,  it  is  known  that  the  re¬ 
duced  characteristic  length  £/Rq  is  a  function  only  of  reduced  concentration  BoC,  where 
C  is  the  concentration,  and  Bo  is  the  second  virial  coefficient.  The  reduced  characteristic 
length  is  the  ratio  of  the  characteristic  length  £,  determined  by  the  angular  dependence 
of  scattered  radiation,  to  the  radius  of  gyration  of  an  isolated  coil,  Rq.  At  low  concen¬ 
tration  £  =  Rq/\/3-  It  is  also  known  that  the  product  of  the  molecular  weight,  M\v,  and 
osmotic  susceptibility  is  also  a  function  only  of  B2C  for  appropriate  solutions.  The 

quantity  can  be  obtained  Uom  the  absolute  scattered  intensity.  Universal  behavior 

is  observed1  over  a  wide  c... -  .utration  range  independent  of  either  molecular  weight  or 
the  solvent  quality  for  polystyrene  both  in  Toluene  (good  solvent)  and  in  MEK  (marginal 
solvent).  Furthermore,  real-space  renormalization  group  calculations2  using  no  adjustable 
parameters  agree  with  these  findings.  The  aim  of  the  present  work  is  to  extend  the  study 
of  the  universality  of  linear  polymers  to  polymers  in  9  solvents  by  varying  the  temperature 
(T)  of  the  polystyrene-cyclohexane  system  from  T  9  to  T  =  9,  at  which  Bo  vanishes.  We 
are  currently  investigating  the  crossover  behavior  from  good  solvent  to  9  solvent  as  tem¬ 
perature  changes,  and  we  hope  to  quantify  the  breakdown  of  universality  as  T  approaches 


lP.  Wiltzius,  H.R.  Haller,  D.S.  Canned  and  D.W.  Schaefer,  Phys.  Rev.  Lett.  51,  1183 
(1983). 

2T.  Ohta,  A.  Nakanishi,  J.  Phys.  A  16,  4155  (1983);  L.  Schaefer,  Macromolecules  17, 
1357  (1984);  T.  Ohta,  Y.  Oono,  Phys.  Lett.  89A,  460  (1982). 

’This  research  is  supported  by  National  Science  Foundation  Grant  No.  DMR-87-14591. 
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"Two-dimensional  reactive  systems:  Rapid  bimolecular 
processes  on  spherical  surfaces" 

by 
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Abstract 
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Crossover  from  Heisenbera  to  Isina  critical  behaviour 
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Shubha  Tewari  and  Sudip  Chakravarty 
Department  of  Physics 
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Los  Angeles,  CA  90024 


Abstract 

It  is  known  from  experiments  done  on  (  A'tjV  1F4),  a  planar  antiferromagnet,  that  there 
exists  a  crossover  from  2D  Heisenberg  to  2D  Istng  critical  behaviour  above  the  3D  Neel 
transition  temperature.  Using  a  nearest  neighbour  Heisenberg  model  with  anisotropy, 
we  do  a  spin  wave  calculation  and  estimate  the  anisotropy  from  a  comparison  with  mea¬ 
surements  of  the  spin  wave  gap.  We  then  map  this  on  to  a  quantum  non-linear  ir-model 
in  2+1  dimensions,  which  describes  the  long  wavelength,  low  temperature  properties 
of  thus  system.  Integrating  out  quantum  fluctuations  we  obtain  an  effective  classical 
model  with  renormalised  parameters.  The  Heisenberg  to  Ising  crossover  temperature 
is  determined.  A  renormalisation  group  calculation  is  done  to  determine  the  crossover 
correlation  length. 
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Curvature  elasticity  of  diblock  copolymer  monolayers 
and  bilavers.  Zhen-Gang  Wang,  Department  of  Chemistry, 
UCLA  —  Curvature  expansion  of  the  free  energy  for  diblock 
copolymers  aggregated  into  monolayers  and  bilayers,  is 
studied.  The  theorv  makes  use  of  the  recent  results  for  grafted 
chains  obtained  by  Milner  et  al.,  and  allows  variation  in  the 
area  per  chain  molecule  with  curvatures.  Elastic  coefficients 
such  as  the  bending  modulus,  the  Gaussian  bending  modulus 
and  the  spontaneous  curvature,  are  calculated  for  a  number  of 
situations  involving  different  conditions  for  the  chains  (melt  or 
swollen).  In  the  case  of  bilayers,  we  show  that,  if  each 
monolayer  has  a  strong  tendency  to  curve  convexly  towards 
the  solvent,  the  flat  geometry  can  become  unstable;  the 
curvature  free  energy  rather  favors  the  formation  of  spherical 
vesicles  over  the  flat  bilayer.  Relevence  of  this  study  to 
microemulsion  using  diblock  copolymer  as  compatilizer,  is 
discussed. 
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A  statistical  mechanical  approach  to  fracture.  Robin 
Blnmberg  Selinger.  Zhen-Gang  Wang.  William  M.  Gelbart. 
Department  of  Chemistry.  UCLA.  —  We  present  a  quasi¬ 
equilibrium  model  of  fracture  of  solid  materials  at  finite 
temperature  and  under  uni-directional  tension.  Our  approach 

is  motivated  by  the  assumption  that  a  solid  under  stress  is 
metastable  thermodynamically.  To  maintain  this  metastable 
state,  tension  must  be  applied  at  a  rate  fast  enough  so  that 
thermodynamic  phase  transitions  such  as  evaporation  and 
melting  do  not  take  place,  and  yet  slow  enough  so  that  local 
molecular  rearrangement  is  possible  for  local  ergodicity  as 
required  by  the  statistical  thermodynamic  approach.  The  solid 
is  modeled  as  an  elastic  bond  network  and  fracture  is 
envisioned  as  the  proliferation  and  aggregation  of  broken 

bonds  under  tension.  We  solve  a  one-dimensional  model 
exactly  and  carry  out  Monte  Carlo  simulations  for  a  two- 

dimensional  network.  Characteristic  properties  such  as  the 

stress-strain  curve  and  its  dependence  on  temperature  are 
obtained.  We  also  discuss  our  results  in  connection  with  the 
Griffith  theory  for  fracture  and  theory  of  nucleation. 
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Abstract 

We  present  an  exactly  solvable  model  of  dissipative  quantum  tunneling  in  a  two- level 
system  coupled  with  a  phonon  bath.  The  model  is  motivated  by  the  problem  of  excimer  for¬ 
mation  in  dimerized  molecular  crystals  such  as  a-perylene.  The  spectral  function  shows  the 
system  in  the  one  dimensional  case  to  be  of  the  "ohmic”  type.  The  decay  behavior  for  this 
dissipative  tunneling  model  may  be  Gaussian,  power-law  and/or  exponential,  depending  on  the 
coupling  strength  and  the  initial  temperature  of  the  lattice. 
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Abstract 

We  study  the  rise  of  an  infinitely  long  bubble  in  a  vertical  cylindrical  tube  full  of  inviscid 
fluid.  We  demonstrate  numerically  that  the  rising  velocity  of  the  bubble  is  determined  by 
the  surface  tension  through  a  solvability  mechanism.  Our  numerical  results  are  in  good 
agreement  with  the  experimental  ones. 


